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A generalization of the Lieb-Schultz-Mattis theorem to higher dimensional spin systems is shown. 
The physical motivation for the result is that such spin systems typically either have long-range 
order, in which case there are gapless modes, or have only short-range correlations, in which case 
there are topological excitations. The result uses a set of loop operators, analogous to those used in 
gauge theories, defined in terms of the spin operators of the theory. We also obtain various cluster 
bounds on expectation values for gapped systems. These bounds are used, under the assumption of 
a gap, to rule out the first case of long-range order, after which we show the existence of a topological 
excitation. Compared to the ground state, the topologically excited state has, up to a small error, 
the same expectation values for all operators acting within any local region, but it has a different 
momentum. 



I. INTRODUCTION 

Lieb, Schultz, and Mattis (LSM) proved in 1961 that 
a one-dimensional periodic chain of length L, with half- 
integer spin per unit cell, has an excitation gap bounded 
by const. /L [1]. This behavior contrasts with the possi- 
bility of a Haldane gap in the integer spin case [2] . 

Despite several attempts [3,4], this theorem has not 
been extended to higher dimensions. The basic difficulty 
in obtaining a higher-dimensional version of this theo- 
rem was pointed out in two insightful papers by Misguich 
and coworkers [5]: if spin correlations are short-ranged, 
the ground state wavefunction should be well described 
by a short-range resonating valence bond (RVB) state 
[6]. The short-range RVB basis decomposes into differ- 
ent topological sectors, depending upon the number of 
dimers crossing a given line through the system. This al- 
lows the construction of a low energy excited state very 
similar to the twisted state of LSM [7]. Instead, if spin 
correlations are long-ranged, such a state will not be low 
energy, but there will exist low energy spin wave excita- 
tions. In contrast to the one-dimensional case, there now 
exist two distinct means of obtaining a low energy exci- 
tation, significantly complicating the proof of any such 
theorem. 

In the present paper, we show a higher-dimensional 
version of the LSM theorem. We consider a d- 
dimensional system of spin-1/2 spins, with finite-range, 
SU(2) invariant Hamiltonian Ti, and with an odd num- 
ber of spins per unit cell on the lattice. Define the total 
number of unit cells in the lattice to be V. Let L be 
the number of unit cells in one particular direction, and 
let L be even; this direction will be referred to as the 
length. Therefore, V is even (if V were odd, there would 
be a trivial spin degeneracy). Let the system be peri- 
odic and translationally invariant in the length direction. 
Let V/L d be bounded by a constant r (this constant r 
is arbitrary, and imposes some bound on the behavior 
of the aspect ratio of the system). Define V/L to be 
the "width" of the system, and let this number be odd. 



Then, we show that if the ground state is unique, the gap 
AE to the first excited state satisfies, 



AE < clog(L)/L, 



(1) 



where the constant c depends on Ti, d, and where the 
result holds for all L greater than some minimum 
where Lq depends on H.,d,r [8]. 

In this paper, we use the term gap to deal specifically 
with the difference between the energy of the first ex- 
cited state and the energy of the ground state. This 
includes two completely distinct physical cases. In the 
case of a one-dimensional system, a spin-1/2 Heisenberg 
chain has a continuous spectrum of excitations above the 
ground state. On the other hand, a Majumdar-Ghosh 
[9] chain has a doubly degenerate ground state with a 
gap to the next excited state. Weak perturbations of the 
Majumdar-Ghosh Hamiltonian can break the exact de- 
generacy between the two lowest states, leaving a system 
with a gap from the ground state to the first excited state 
which is exponentially small in the system size, and then 
a gap from the first excited state to the next excited state 
which is non- vanishing even in the limit of large system 
sizes. We consider both of these cases as systems in which 
the gap AE is vanishing in the limit of large system size. 
Although they are one-dimensional systems, these two 
cases closely match the two possibilities mentioned above 
for higher-dimensional systems. The first case involves a 
system with a continuous spectrum as it has algebraically 
decaying spin correlations. In the second case, the first 
excited state is very close to the twisted state of LSM. 

The physical idea behind the proof of Eq. (1) is closely 
related to the two possibilities considered above for the 
absence of a gap. In the event of long-range order, or 
algebraic long-range order, one expects that there is no 
gap. Conversely, if there is a gap, one expects that there 
is no long-range order. This is the first statement we 
prove: we assume that the system has a gap AE and 
show, in section III, that connected expectation values 
decay exponentially in the spacing between them. Then, 
to prove Eq. (1) we first assume that Eq. (1) is violated, 
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proceeding by contradiction. We assume the existence 
of such a gap violating Eq. (1) in sections IV, V, and 
VI. We then use the existence of the gap AE and the 
exponential decay of correlation functions to show an in- 
sensitivity of the system to boundary conditions in sec- 
tion IV. Then, this insensitivity is used to construct a 
low-energy, twisted state in section V. The construction 
of this twisted state will to some extent follow the topo- 
logical attempt [4] at proving the LSM theorem in higher 
dimensions, with some important differences outlined be- 
low. In section VI we will show that the twisted state has 
a different momentum than the ground state. It is here 
that the odd width of the system becomes essential. De- 
spite the different momentum compared to the ground 
state, the twisted state has, up to a small error, the 
same expectation values for all operators acting within 
any local region. Thus, we may refer to this state as a 
topologically excited state. Finally, in an appendix, we 
briefly consider a version of the result showing exponen- 
tial decay of correlation functions for systems governed 
by certain Markov processes, rather than quantum sys- 
tems. 

Since we will constantly deal with operator equations 
of motion, we introduce a set of "loop operators" as a ba- 
sic technique. The loop operators, which will be suitably 
defined products of spin operators, can be naturally in- 
terpreted as a product of gauge fields around a loop [10]. 
However, the use of these operators avoids the uncon- 
trolled approximations associated with the U(l) and Z-i 
gauge theory techniques [11]. The introduction of these 
operators is not necessary to the main development, but 
provides a useful notation. 



operator tr(ijk...m) permutes the spins around the sites 

Given an operator 0(t), the operator obeys the equa- 
tion of motion dtO(t) = —i[0,TC\. Consider, for example 
a term tr(ij') in TL. We have [i,tr(r?')] = ij — ji. As 
an illustration, let us give the full Greek indices on this 
commutator: we have [i» v , {i a Pf a )] = i^ a j av - j^ a i av 
where summation over repeated indices a, (3 is assumed. 

Introduce coordinates (x, y) to specify sites i, where 
x labels the unit cells along the direction of length and 
is defined up to integer multiples of L. The coordinate 
y labels the unit cells along the other lattice directions, 
as well as labeling the particular spin within the unit 
cell. Given two sites i and j on the lattice, we define 
the distance between them, written \i — j\, as the mini- 
mum number of moves by lattice vectors needed to move 
from the unit cell containing i to that containing j. On 
a square lattice, for example, this is the Manhattan dis- 
tance. Then, let R denote the range of Ti, the furthest 
distance between two sites in any term in H. If all dis- 
tances in a product of loops are less than L/2, we can 
define a winding number of the given product around 
the lattice in the length direction. If all distances in 
the product remain less than L/2 — R, then the dynam- 
ics dtO — —i[0,H] does not connect sectors with differ- 
ent winding numbers. We will make use of the coordi- 
nates later by sometimes writing loop operators tr(«j...) 
in the form tr((a;i, yi)(x2, 2/2)...), where i has coordinates 
(x\,yi), j has coordinates (£2,2/2), and so on. 



III. LOCALITY 



II. LOOP ALGEBRA 



We define operators = 



, £a W. where 
Sf are the spin operators at site i and a a are the Pauli 
matrices, a — x,y, z. Thus, is the two-by-two matrix 
of spin operators 



s? + , 



Sf s?- l s? 

1 QZ 

2 S i 



(2) 



We consider operators of the form V 11 'j vp k pa ...m afl , 
which we refer to as loop operators, where a summa- 
tion over repeated indices is implied. Later we will often 
suppress the indices fi, v, writing i rather than V 1 " to save 
space. Thus, we will write the loop operator mentioned 
above in the form tr(ijk...m), where the trace tr refers 
to a trace over the Greek indices /1, v, .... Below we also 
use a trace Tr; this trace Tr refers to a trace of quantum 
operators, summing over all states in the Hilbert space of 
the system. Using the rule i 9V i pa = 6 ,yp i fJ '' 7 , it is always 
possible to reduce a given product of traces to a new 
product such that each site appears only once. Then, an 



We consider ground state expectation values of oper- 
ators Oi,0 2 ---, written (Oi(ii)0 2 (i2)---)- The expecta- 
tion values are not time ordered: the ordering of oper- 
ators is as written. For a system with a unique ground 
state and an energy gap AE, on physical grounds one 
expects that connected correlation functions, defined as 
(A{0)B(0)) c = (A(O)S(O)) - (A(0))(B(0)}, decay expo- 
nentially in distance (without loss of generality, we will 
assume (A) = (B) = through the rest of this section). 
The proof of this locality bound will be done in this sec- 
tion. Wc will do this in two steps: first, we consider com- 
mutators of the form [A(t),B(0)}, where A(0) and B(0) 
are separated in space. We bound the operator norm 
[13] of the commutator for sufficiently small t, and thus 
bound its expectation value, in Eq. (4) below. The proof 
in this subsection will just be sketched; a more rigorous 
derivation is due to Lieb and Robinson [14]. This result 
provides a bound on the velocity of the system, as will 
be seen below. Then, in the next subsection, from this 
bound on the expectation value of the commutator and 
the existence of a gap, we use a spectral representation 
of the commutator to bound the connected correlation 
functions, thus obtaining the desired locality bound on 
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the expectation value, Eq. (13). Finally, we close the 
section by giving a similar locality bound for operators 
separated in time. 



A. Finite Velocity 

We define the distance between two operators 0\ , O2 
to be I if the minimum distance between any pair of sites, 
where Si appears in 0\ and Sj appears in 2l is I- 
A(0), -B(O) are sums of products of spin or loop operators, 
which wc suppose to be distance I apart. 

We start with some notation. The Hamiltonian, H, 
can be written as a sum of terms H = ^ Hi, such that 
each Hi only contains spins operators on sites j with 
\i — j\ < R- Let Na denote the number of sites appearing 
in A(0), and Nb denote the number of sites appearing in 
B(0). Let J denote the maximum, over sites i, of \\Tii\\. 

We now bound the operator norm of the commutator 
[A(t), -B(O)] for short times. On short time scales, one 
expects that A(t),B(0) are still separated in space, up to 
small correction terms, as we now show. Consider first 
||[A(f), "Hi]\\, and study the change in this quantity as a 
function of time: 



* E 

\i-j\<2R 

- E 

\i-j\<2R 



\\[A(t)-idt[A(t),H],Hi]\\-\\[A(t),Hj 
dt 

\\[A(t),Hi + idt[Hi,H]]\\ - \\[A(t),Hj 
' 1 dt 

A\[A(t),H i + idt[H i ,H j ]]\\-\\lMt),'Hi 
dt 

WlA^-idtlA^^lHM-WlMt)^ 
dt 

<2J lUiWj 

\i-j\<2R 



(3) 



Here, we work to linear order in dt. While the operator 
[A(t),Hj] is differcntiable, its operator norm need not be. 
Thus, all equations here are correct when we take the lim 
sup as dt — > 0. 

The first equality in Eq. (3) is obtained by moving 
the time derivative from A(t) to Hi as follows: for 
any operator P, to linear order in dt we have ||P|| = 
\\P + idt[P,H]\\. Set P = [A(t)-idt[A(t),H],Hi\. Then, 
to linear order in dt, \\P\\ = \ \P + idt[P, H}\\ = \\[A(t) - 
idt[A(t),H],Hi] + idt[[A(t),Hi],H}\\ = \\[A{t),Hi + 
idt[Hi,H]}\\. 

The inequality is obtained because [Hi,Hj] = for 
\i — j\ > 2R. The next equality is obtained by moving 
the time derivative back to A(t), using now the equality 
= idt[P, Hj]\\. The final inequality results from 
the bound \\Hi\\ < J. 

Now, let S denote the maximum number of sites j 
within distance R of any site i. Eq. (3) gives a set of 
differential equations which bound the operator norm 



of various commutators; we have also the initial con- 
ditions that I |[A(0),7<j] 1 1 vanishes for sites j which are 
further than distance R from any site in A(0), while 
||[A(0),Wj]|| < 2J||A|| for all other sites. The number 
of sites within distance R of A(0) is bounded by NaS. 

To bound ||[A(i), let us then consider the fol- 

lowing set of differential equations: for t > 0, we 
take d t Gi = 2JJ2\i-j\<2R.Gj and for t < we take 
d t Gi — —2JJ2u_ji <2R Gj, with initial conditions Gj — 

for sites j which are further than distance R from 
any site in A(0), and Gj = 2J||A|| for all other sites. 
Then, comparing these equations to Eq. (3), we see that 
1 1 [A(t), Hi] 1 1 < Gi(t). This set of linear equations for 
d can be solved for any given lattice. However, we 
are simply interested in an upper bound on d. Let 
us define G k to be the maximum of Gi over all sites 

1 which are at a distance greater than (2k — 1)R from 
all sites in A(0). Then, we have d t G k < 2JSG k ~ 1 
for k > 0, and d t G k < 2JSG k for k = 0, with ini- 
tial conditions G° = 2J||A|| and G k = for k > 1. 
Thus, G°(t) < 2J\\A\\e 2Jst , G^t) < f* dt'2JSG°(t') < 
J *dt'(2J5)(2J||A||)e 2J5t = 2J||A||(2J^)e 2J5t , and 
G k (t) < 2J\\A\\{2JSt) k e IJSt /k\. The last set of inequal- 
ities follows inductively: G k (t) < f* dt , (2JS)G k - 1 (t) < 
f Q dt'{2JS){2J\\A\\){2JSt') k - 1 /{k - l)!e 2JSt 
2J||A||(2JS't) fc e 2J5 7fc!- From these inequalities, wc find, 
for a site j which is at a distance greater than (2k — 1)R 
from A(0), that \\[A(t), Hj}\\ < 2jp|||2JS't| fc e 2Js l*l/fc!. 

Finally, consider <9 t ||[A(i), B(0)]||. Using a simi- 
lar sequence of inequalities to Eq. (3), wc find that 
d t 1 1 [A(t), B(0)} 1 1 < 2| \B\ I £ . 1 1 [A(t) ,Hj]\\, where the sum 
over j extends over sites j which are within distance R 
of some sites in B(0). There are at most NbS such 
sites, and each of them has ||[A(t), Hj]\\ < G fe_1 (t), 
where k = 1/2R. Here, we take k to be the ceiling 
of 1/2R, the smallest integer greater than or equal to 
1/2R; we obtain this value of k since each such site is at 
least a distance I — R from all sites in A(0), and so we 
need I- R > [2(k - 1) - 1]R. Then, || [A(t), B(0)}\\ < 
2N B \\A\ 1 1 \B\ 1 12 JSt\ l ? 2R e 2JS M /(l/2R)\. 

Define f(t) = (\A(t), 5(0)]). Since f(t) < 
||[A(t),B(0)]||, 



/(*) < 



2jVB||A||||g|||2J^|'/ 2fi e 2Js l*l 
(l/2R)\ 
= N B \\A\\\\B\\g(t,l). 



(4) 



For t = Cil, the large I behavior of g(t,l) <~ 
exp[(l/R)(2JS Cl R + 1/2 + (l/2)log(4JS'cii?))]. If we 
choose a sufficiently small Ci, then g(c\l,l) decays ex- 
ponentially in / for large I. Numerically, we find that 
the zero of 2JS Cl R + 1/2 + (l/2)log(4JS'cii?))], is at 
ci w 0.139232/(2JS' J R). Any ci smaller than this value 
(for example, c\ — 0A/(2JSR) will work) will cause 
g(cil,l) to be exponentially decaying for large I. The 



3 



velocity at which correlations spread in the system is of 
order c^ 1 . 

B. Spectral Decomposition 

Now, we use a spectral decomposition of /(f) to re- 
late /(f) to the desired correlation function, (^4(0)5(0)). 
Without loss of generality, let us set the ground state en- 
ergy, £o, to 0. The spectral decomposition of /(f) gives 

/(f) = J2(e- iEit A 0i B i0 - e iEit B 0i A i0 ), (5) 

i 

where Aio is the matrix element of operator A between 
the ground state and the eigenstate i, with energy 
Ei > AE above the ground state energy, and similarly 
for the other A^, B^, Boi. There are no terms in Eq. (5) 
involving A 00 , B 00 since we have assumed (A) = (B) = 0. 

Let us define a function ,f + {t) = J2 i e~ Eit AoiBio, 
which thus contains only the negative frequency (posi- 
tive energy) terms in /(f). The significance of f + (t) is 
that f + (t) = (A(t)B(0)), so that the positive energy part 
of /(f) contains the desired correlation function in it. In 
this subsection, we combine the bound (4) on /(f) with 
the existence of a gap to bound .f + (0). 

Define /(f) = f(t)e- t2AE ^ (2c i\ with q to be cho- 
sen later. We have two bounds on /(f). First, we 
have the bound (4) on /(f) which gives us the bound 
/(f) < N B \\A\\\\B\\g(t,l). We also have have 

/(f)<2P||||B||e-* 2Aij2 /^. (6) 

We will use the first of these bounds for times |f| < c\l, 
and the second for long times |f | > c\l. Finally, we define 
f + {t) to contain only the negative frequency terms in 
/(*)■ 

Now, the desired expectation value (A(0)B(0)) = 
/+(0). To bound /+(0), we first bound /+(0), and 

then bound /(0) + - /+(0). To bound /+(0), we use 
the bounds on / and an integral representation of the 
positive energy part [15,16]: 

1 r°° 

l/+(0)l = 2^' J J t f( t )/(-^ + ^)\ 

= ^(\[ Kt)/{-it + e)\ + \ f f(t)/(-it + e)\) 
llT v J\t\< Cl l J\t\> Cl l ' 

<^||A||||i?||(2^(c 1 M) + 2^ e -^ 2As2 /^). 

In Eq. (7), to bound the integral over |f| < C\l, 
we used |/(f)| < |/(f)| < JV B ||A||||B||(/(t,0 < 
Ar B ||A||||B||[|f|/(ci0]5(ci?,0- To derive this inequality 
we have assumed that I > so that taking the ceiling of 
1/2R above gives a k > 1. Then, 



. / df/(f)/Hf + e)| < (8) 
J\t\<al 

N B \\A\\\\B\\ [ dtg(c 1 l,l)/(c 1 l) = 
J\t\<al 

2N B \\A\\\\B\\g( Cl l,l). 

To bound the integral over |f| > c\l in Eq. (7) 
we have used Eq. (6) to show | Jui >Cl ; f(t)/(— it + 
e)\ < 2||A||||B||/_ oo oo dfexp[-f 2 A J B 2 /(29)]/(ci0 < 
2\\A\\\\B\\j^eM~clPAEy(2 q )}. 

To bound |/(0) - /(0) + |, we start with the definition 
of /. Expressed as a convolution in Fourier space [17] 
this is: 

/(a/) = (v/2^/A£) J dLjf(Lu)e- q ^'^^ 2AE2 K (9) 

Now is where the existence of an energy gap be- 
comes essential. For motivation, let us first pictorially 
(see Fig. 1) describe how the gap enables us to bound 
|/ + (0) — / + (0)| and then present it more mathemati- 
cally. By definition /(0) = /(0); this follows in Fourier 
space from J M dw/(w) = J" dw/(w). The convolution 
(9) means that a given Fourier component in / which is, 
for example, negative frequency, will produce both pos- 
itive and negative frequency Fourier components in /. 
So, consider a (5-function spike in f(u>), produced by an 
intermediate state i with energy Ei = —ui > 0. This 
produces a Gaussian in /(w), as shown. The integral 
over all ui of the Gaussian is the same as the integral 
of the J-function; however, the shaded portion of the 
curve has to > 0. Since /+(0) = (27T)" 1 j"^ duf(u) 
and / + (0) = (2n)^ 1 j"^ duf(u)), we find a difference 

between / + (0) and / + (0) equal to the integral of the 
shaded portion of the curve. At ui — the height of the 
Gaussian is reduced by a factor er qu /( 2AE ). However, 
since > AE, this factor is bounded by e~ q / 2 . 

Now, let us do the calculation more directly: / + (0) = 
(27T)" 1 du>f{u>) du'iV^/AE) exp[- g (o; - 
co') 2 /(2AE 2 )}, while /+(0) = (2^ J° x dcuf(uj). Then 

(7/+(o) - /+(o) = (27T)- 1 r d w /(c)[e 9 (-c) - e(- w )]. 

J — oo 

(10) 

Here 8(w) is a step function: 0(w) = 1 for uj > and 
O(w) = for lo < 0. We have defined 

/>oo 

0,(w)= / dJ{y/2^q/^E)exp[-q(u- J) 2 /{2AE 2 )). 
Jo 

(11) 

Since the system has a gap, the integral in Eq. (10) van- 
ishes for \uj\ < AE. However, for \u)\ > AE, we have 
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\Q g (-uj) - Q{-lo)\ < e- q ' 2 . Thus, Eq. (10) is bounded 
by 

/oo 
dc|.fH|<2||A||||B|| e ^/ 2 . (12) 
-oo 



Thus, combining Eqs. 
|/+(0) - /+(0)| + |/+(0)| < 

2 ^||I e -^' 2 AB 2 /(2 9 )) + 2P||||B|| e -9/ 2 . 

choose q = cilAE to get 



(7,12), |/+(0)| < 
±\\A\\\\B\\(2N B g( Cl l,l) + 



We finally 



|(A(0)B(0)) c | < —2N B \\A\\\\B\\g( Cl l,l) + 



(13) 



2(1 



me - Cl lAE/2 



giving the desired bound. The first term in Eq. (13) 
decays as exp[— 0(l/R)], while the second term decays 
as exp[-ci/A£;/2] = cxp{—0[AEl/(JSR)]}; here, by 
0(l/R), we mean some quantity of order l/R. 

In what follows in the next three sections, the first term 
in Eq. (13) will be negligible: we will be considering oper- 
ators separated by a distance I which is of order L, so that 
the first term in Eq. (13) will lead to only exponentially 
small (in L) contributions to the correlation functions. 
The second term will be more important: since we will 
consider gaps AE oc log(L)/L, the second term will lead 
to terms which are suppressed only by powers of L when 
considering correlation functions of operators separated 
by a distance of order L. 



C. Operators at Different Times 

It is possible to extend the result Eq. (13) to cor- 
relation functions (A(O)B(iT)), with r real and r > 
0. Then, in Eq. (7), we must evaluate |/ + (— ir)\, so 
that the denominator (—it + e) is replaced by — it + 
t. In this case, we are still able to find just as 
tight a bound on |/ + (— ir)\ as we previously found 



< 27 



_ e ^?i 2 A£ 2 /(2?) N | 

J A£ci( c I' 



B\\(2N B g( Cl l,l) + 



for |/+(0)|: \f+(-ir 

Of course, for r > qj 'AE there is the trivial 
bound \(A{0)B(it)) c \ < ||A||||B||e^ AE < ||^||||B||e-«. 
For |t| < q/AE, we claim that \f+(-ir) - 
cx V [+T 2 AE 2 /(2q)]f+(-iT)\ < 2P|||| J B||e-9/ 2 . To show 
this, \f+(-iT)-cx V [+T 2 &E 2 /(2q)]f+(-iT)\ = 



(27T)- 1 / dwf(u)([ AJ{y/2^/AE) 

J — oo J — OO 

exp[u/r] ex-p[-q(io - lo') 2 / (2AE 2 )] - 
exp[+T 2 AS 2 /(2g)]6(-w) cx P [ut]) . 

The portion of the integral with to < is equal to 



(14) 



<*)-/_ 



-AE 



Aujfiuj) 



(15) 



exp[u/r] exp[-q(io - lo') 2 /{2AE 2 )], 



where we have used the gap AE and the relation 

&Lo'(*fh^/AE) exp[u/r] cx V [-q(Lo - u') 2 /(2AE 2 )] 
= exp[wr] exp[+r 2 AS 2 /(2(7)]. Then, for r < q/AE, the 
integral (15) with lo < is bounded in absolute value by 
(2tt)^ 1 f_ dto\f(io)\e~ q / 2 . We can similarly bound the 
portion of the integral with lo > 0, giving the desired 
result. 

With the given q = c\lAE the above bounds show that 
for r < ciZ, \(A(0)B(iT)) c \ < 



-T 2 AB/(2cii)r 



2(1 



l 2tt 
1 



^/2ir Cl lAE 



2N B \\A\\\\B\\g(c 1 l,l) + 
)||A||||B||e- ClAra ' /2 ]. 



(16) 



IV. TWISTED BOUNDARY CONDITIONS 

In this section we derive some results on the sensitivity 
to boundary conditions, as a step towards the the main 
result, Eq. (1). To derive a contradiction later, we will 
assume throughout this and the next two sections that 
there is a gap AE that violates Eq. (1), with an appro- 
priately chosen c. In the first subsection, we review the 
twist of boundary conditions and the topological attempt 
at proving the LSM theorem. In the second subsection, 
we show the specific results on the sensitivity to bound- 
ary conditions. 



A. Topological Argument 

Here we will define a new twisted Hamilto- 
nian, making use of the coordinates, (x,y), in- 
troduced previously for lattice sites i. To de- 
fine the new twisted Hamiltonian, Tig^i, replace all 
loop operators tr((xi, yi)(x2, y2)(x 3 , 2/3)...) in Ji with 
tr((x 1 ,y 1 )R(x 1 ,X2){x2,y2)R{x2,x 3 )(x 3 ,y 3 )...). Here, 
the twist operator R tJ - u (xi,x 2 ) = cxp[±i^a^] if the 
shortest lattice path between xi,X2 crosses from x = 
to x = 1, where the sign is positive if the path crosses in 
the direction of increasing x and negative if it crosses in 
the opposite direction. Here, R fl ' y (xi,X2) is a two-by-two 
matrix of numbers, rather than of operators, 



BP" = 



exp[±i|] 
exp[=Fi|] 



(17) 



Alternately, if the shortest lattice path between x\,X2 



crosses from L/2 toi/2 + 1, R^(x 1 ,x 2 ) = exp[±i^-a^}. 
Otherwise, R(xi,x 2 ) = 5^ v . In Fig. (2), we show the 
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coordinate system using x, y and show where the two 
boundary condition twists are inserted. 

Let us see what the effect of this twist is in terms of spin 
operators. Consider two sites, i, j. Suppose the Hamil- 
tonan Ho,o has a term such as tr(ij) = 2(SfSj + SfS* + 
SfSj ) + 1/2. Then, let us suppose i has x = while j 
has x = 1. Then, He,e' has a term tr(ii?(0, l)jR(l, 0)). 
In terms of spin operators, this is equal to 2[S?Sj + 
cos(6){S?S? + SfS]) + sm(6){S?S? - SfS*)} + 1/2. In 
the untwisted Hamiltonian, we coupled the dot product 
of the two spin vectors, £>i,Sj; in the twisted Hamilto- 
nian, we couple them after rotating one by an angle 9 
about the z-axis. A good discussion of twists can be 
found in [5]. 

We have considered two different twist angles, 9,9'. 
The spectrum of Tlg y g> depends only on the combination 
6 + 6' . Further, from any given eigenfunction "0(0, 9') of 
Hg.e>, one can find an eigenfunction ip(9 — 89,9' + 59) 
of Hesefi'+se by - 60, 9' + 69) = ft. e ms ^(6, 9'), 
where the product extends over all sites j = (x, y) with 
< x < L/2. 

Given that the spectrum depends only on the combi- 
nation 9 + 9' , the reader may wonder what the reason 
is for introducing two twist angles, rather than just one 
angle. In fact, the second angle is a useful trick, intro- 
duced for the following reason: we have previously shown 
that the existence of a gap causes correlation functions to 
decay exponentially in the separation of the two opera- 
tors. However, physically, one expects that the existence 
of a gap will also imply some insensitivity of the system 
to boundary conditions, enabling us to bound, for ex- 
ample, the second derivative of the ground state energy 
with respect to 9. What we will do in the next subsection 
is show this insensitivity by using the fact that the spec- 
trum depends only on 9 + 9' to convert the second deriva- 
tive (dg) of the ground state energy into a mixed partial 
derivative (dgdgi) of the ground state energy, and by then 
evaluating that mixed partial derivative as a correlation 
function, using the exponential decay of correlation func- 
tions. This will be stated more precisely at the start of 
the next subsection; we mention it here for motivation. 

The eigenvalues are of Hg,gi are invariant under 9 + 
9' — > 9 + 9' + 2tt, while the wavefunctions are invariant 
under 6 + 6' -> 6 + 0' + 4ir [4,5]. To motivate the results 
in this section, we recall the basic idea of the topolog- 
ical attempt [4] at proving the LSM. The idea is that 
if there is a gap at 6 + 9' = 0, and if the gap remains 
open for all 9 + 9', then under an adiabatic change in the 
angle 9 with 9' fixed at zero, the ground state at 9 = 
evolves into the ground state at 9 — 2n. At 9 = 2ir, the 
Hamiltonian is returned to the original Hamiltonian, but, 
for a system of odd width, the ground state expectation 
value of the translation operator changes sign, as will be 
discussed in more detail below. This leads to a contradic- 
tion: from the ground state with given expectation value 



of the translation operator, we construct another ground 
state with the opposite expectation value. The require- 
ment that the topological attempt requires the gap to 
remain open for all 9 was pointed out in [5]. 

What the topological argument actually succeeds in 
showing is that the gap must close at some value of 9. 
However, in order to use this argument to obtain any 
bound on the magnitude of the gap at 9 = 0, we would 
have to show that a sufficiently large gap at 9 = would 
prevent the gap from closing for all 9; that would then 
lead to a contradiction, enabling us to bound the gap at 
9 = 0. What we will see is that we can partially show 
this: for sufficiently large c in Eq. (1), we can show to 
second order in 9 (or indeed, to any finite order) a bound 
on the change in ground state energy with respect to 9. 
However, we will be unable to show that the gap remains 
open for all 9 because to bound the change in ground 
state energy for higher orders in 9 requires progressively 
increasing the constant c in Eq. (1), and it is not possible 
to show the result to all orders. Thus, the topological 
attempt will ultimately fail, and we will give a physical 
example of how this can happen. In the next section (V), 
we will give a successful argument. 

B. Boundary Condition Sensitivity 

We now show an insensitivity of the ground state en- 
ergy, Eq{6,6'), to second order [18] in the twist angle, 
9 + 9'. At 9 = 9' = 0, d 9 E (9,9') = {d e H e ,e<) = 0. In- 
deed, taking any odd number of derivatives of Eq(9,9') 
leads to a vanishing quantity [19]. To second order in 
9,9', we write a power series: E (9,9') = E (0,0) + 
a9 2 /2 + a9' 2 /2 + b99' , where a = d^E = d 2 ,E and 
b = Ooo'Eq. We will show that, for any given negative 
power of L, we can find a constant c such if Eq. (1) is vi- 
olated for that c, then a is bounded by an 7i-dependent 
constant times the given negative power of L. We do 
this by calculating b as a correlation function, and then 
showing that b = a. 

Recall linear perturbation theory: suppose a Hamilto- 
nian 7i is changed by some 5TL. For a non-degenerate 
state, \ip), with eigenvalue E, the change \8ip) in \ip) is 
given to linear order in STi. by \5ip) = (E — H)~ 1 5H\ip). 
Since the ground state is the lowest energy state, all 
other states have energies greater than it. Thus, we 
can write the change in the ground state to linear order 
as \SJ ) = -E a ^ir dTeiEo ~ Ea)T \^)(^\ Sn ^o) = 
— L dTdH.(iT)\ipo) , where |^o) is the ground state wave- 
function, \ip a ) are a complete set of intermediate states, 
and where SH(it) = cxp[— Ht]SH cxp[Hr} is the change 
in the Hamiltonian operator, taken at imaginary time it. 
Here we have set E o (0, 0) = without loss of generality. 

Specializing to the case of SH = dgHefi 1 and writing 
the change in ip in terms of the 9, ^'-dependent ground 
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state density matrix p°{6,6') = \ip (9,6')}(ipo(6,9')\ we 
have: 

/>oo />0 

dg P ° = - &Td e H e ,e>(iT)p - / dr p°d e He,e>{iT), 

JO J-oo 

d >p° = - dTd >He,e>(iT)p° - / dr p°d e ,H e fi>{iT). 

JO J-oo 

Note that since (<57i) vanishes in this case, we do not need 
to worry about matrix elements of SH from the ground 
state to the ground state. 

Now, we can use the change in the density matrix to 
compute b by b = Tr[(d e >H e ,e>){deP )]- So, 

b = -J o dT ((deH e ,e>{0)de>He,e>{iT)) (19) 

-(do>He, >(-iT)d e Ho,o>(0))), 

where the derivatives are evaluated at 6 = 0' = 0. The 
derivative dgHi is non- vanishing only for sites i which are 
within distance R of x = 0; there are at most SV/L such 
sites. For each i, \\d e Hi\\ < JS, so \\d e H e ,e>\\ < JS 2 V/L. 
We use two bounds for the given correlation functions in 
Eq. (19). First, each correlation function is bounded by 
( JS 2 V/L) 2 e~ TAE . Second, we can use Eq. (16) to bound 
each correlation function by 



2S(JS 2 V/L) 2 {1 + l/y/Trc 1 LAE)cx.p[-c 1 AE(L/2)/2], 

(20) 

where we neglect the term in g(c\l, I) in Eq. (16) as it 
leads to a correction which is exponentially decaying in 
L, not in c\AEL, and thus is negligible in what follows. 
Also, we have used I = L/2, ignoring the slight error, 
that in fact I > L/2 — R. Finally, we have used Nb < S 
in Eq. (20). 

Using these two bounds on the correlation function, we 
arrive at 

oo 

|6| < 2( JS 2 V/L) 2 j drmin{exp[-rA£;],a;}, (21) 



where x = 25(1 + 1 /s/ira LAE) exp [- a AE(L /2) /2] . 
Thus, |6| < 2(JS 2 V/L) 2 (x/AE){1 +log{x)). The num- 
ber of sites, V, is bounded by rL d , while for AE greater 
than c\n{L)/L, x is bounded by a c-dependent negative 
power of L. Therefore, we can bound |6| by any desired 
negative power of L by choosing c sufficiently large. 

However, E (9,-9) = E o {0,0), so b = a. Thus, we 
have also bounded |a| by the same negative power of L. 
Therefore, at 6 = 6 1 = 0, we find that \d$E (6,9')\ is 
bounded by a negative power of L. This shows some in- 
sensitivity of the ground state energy to boundary con- 
ditions. This realizes the physical idea [5] that a spin 
liquid state is defined by the lack of response to a twist 
in boundary conditions to second order in 9 [20]. 



At fourth order in 9, we must evaluate a correlation 
function of jour operators, each of order JS 2 V/L; to 
bound these correlation functions requires a larger c. 
(lE)ach higher order in 9, 9' requires an even larger c, so 
that it is not possible to bound the change in ground 
state energy for arbitrary 9 + 9' . Therefore, the topologi- 
cal attempt [4] to establish the LSM result fails. Indeed, 
a gap at 9 + 9' = must close for 9 + 9' ^ [5]. 

It is worth giving a specific physical example of this 
possibility, as the topological argument does show that 
the gap must close for some 9 + 9' . In many physical 
examples of spin liquids, the closing of the gap arises 
because a state which is at some very low energy, of order 
JL^ 1 or less, above the ground state at 9 + 9' = crosses 
the ground state energy at a finite 9 + 9'. For example, if 
the Majumdar-Ghosh Hamiltonian is slightly perturbed, 
there is a state at an exponentially small energy above the 
ground state which crosses the ground state at 6+6' = it. 

However, it is also possible for a state which is at some 
energy JL° to cross the ground state: consider a system 
with two competing phases, one of which is a spin liq- 
uid phase while the other is a spin ordered phase with a 
spiral order. The spiral order is chosen so that the spin 
ordered phase can be frustrated at 9 + 9' = 0, and the 
spin liquid is the ground state there. At some 9 + 9' ^ 0, 
however, the spiral phase can take over as the ground 
state. This taking over as the ground state can happen 
either via a level crossing (if the two states have different 
symmetry, for example, or if the spin ordered phase has a 
non- vanishing net spin) , or via an avoided crossing. This 
provides a specific example of a system in which a state 
or phase which is at an energy of order JL° at 9 + 9' = 
becomes the ground state at some non- vanishing 9 + 9' . 

The solution to this problem is simple: it is not nec- 
essary to show that there is a gap for all twist angles. 
Instead, we start with the ground state at vanishing 
twist and continuously evolve this state, obtaining a state 
for any twist angle which is an approximate cigenstate 
of the twisted Hamiltonian, not necessarily the ground 
state. This approximate cigenstate will be explicitly con- 
structed in the next section, while in the section after 
that we demonstrate that at a twist of 2ir the expecta- 
tion value of the translation operator has changed sign in 
the new state compared to the ground state. Thus, this 
gives a new low energy state, different from the ground 
state. 



V. TWISTING THE GROUND STATE 
A. Constructing the Twisted State 

Let p(6, 0) be a ^-dependent density matrix that we 
construct below. Divide the system into two overlapping 
halves: half (1) contains sites with x = 3L/4 — R, 3L/4 — 
R + 1, L - 1,0, 1, ...L/4 + i?, while half (2) contains 



7 



sites with x = L/A - R, L/A - R + 1, 3L/4 + R. That 
is, half (1) contains all sites from x = 3L/A — R up to 
x = L — 1, as well as all sites from x — up to L/A + R, 
while half (2) contains all sites from x = L/A — i? up to 
x = 3L/A + R. The halves are shown as shaded regions 
in Fig. (3). 

The reason we choose two overlapping halves is that we 
will be considering density matrices which involve only 
sites within a given half. These matrices will be defined 
by tracing over sites outside the given half. Then, to 
evaluate the expectation value of the energy of the sys- 
tem, we will be able to evaluate the expectation value as 
a sum of operators which lie completely within one or the 
other half. That is, by making the two halves overlap, 
we will deal with the question of the "seam" where the 
two halves join. This is mentioned here as motivation 
and will be done in more detail below. 

Define pi(0,O) = Tr 2 [p(6», 0)], where Tr 2 denotes a 
trace over all sites not in half (1), and define p2($,0) = 
Tri[p(0, 0)], the trace over sites not in half (2). Sim- 
ilarly, define p?(0,0') = Tr 2 [p°(0, 9% and p° 2 (6,9') = 
Tri[p°(#, #')]. We will assume throughout this section 
that there is a gap violating Eq. (f). Then, for suffi- 
ciently large c, we will construct pi such that 



p 1 (e,0)-E 1 (6)=(j[e i ^)p 1 (0,0)(j[ t 



(22) 



P°i(0,-0), 



where the products extend over all sites j — (x, y) with 
< x < L/A and where E\{9) is an error term such 
that the trace norm [21] |-Ei(#)| is bounded by a con- 
stant times a negative power of L for all < 9 < 2n. 
The particular negative power of L can be determined 
by choosing the constant c in Eq. (1). As a useful termi- 
nology, we will refer to a quantity as "small" if, for any 
desired negative power of L, we can find sufficiently large 
c or sufficiently large q (introduced below), such that the 
given quantity is bounded by a constant times the given 
negative power of L for all L. Thus, we wish \E\ (9) | to be 
small. Note that, given this definition of small, if a small 
quantity is multiplied by any fixed power of L, the result 
is a small quantity. Sometimes, we will indicate that a 
quantity is made small by choosing c or by choosing q, 
to specify which of the two needs to be made sufficiently 
large. 

In differential form, we require 

fcpi(M) = £i[SM(0,-0)]+ei(0), (23) 



where e\ = dgEi. We will show that the upper [22] 
derivative Vg\Ei(9)\ is small, from which Eq. (22) will 
follow. We will also require p2(9,0) = /^(O, 0), up to 
a similarly bounded error term i?2(#), and dgpi(9,0) = 
e2(9), with a similarly bounded Vg\E2{9)\. 



The physical motivation behind Eq. (22) is to construct 
a state for the Hamiltonian Hg t o that has an energy close 
to E o (0,0). The twist 9 is along a line that lies com- 
pletely within half (1) while 9' is along a line that lies 
completely within half (2). Within half (1), the Hamilto- 
nians Tig^o and Tig-g are equal, so we construct a density 
matrix such that within half (1) the given density ma- 
trix is close to the ground state density matrix of TCe t -$. 
Then, the expectation of any operator O which lies com- 
pletely within half (1) for the density matrix p(9,0) will 
be within ||0|||i±a(#)| of the expectation value of that 
operator for the density matrix p(9 1 —9). On the other 
hand, within half (2), the Hamiltonians Tig_o and 7io,o 
are equal, so we also require that within half (2) the den- 
sity matrix be close to the ground state density matrix 
of H ,o- 

Then, the expectation value of the energy in the state 
defined by p(9, 0) is equal to Tr[p(9, 0)Hgfi}- Once we 
have shown that both Eq. (22) and the bound on |£ l 2(^)| 
are satisfied, it will follow that this expectation value will 
be within an amount ||W0 ; o||max(|.Ei|, I-E2I) 01 ^o(0, 0), 
since the Hamiltonian TL can be written as a sum of oper- 
ators which are entirely within half (1) or entirely within 
half (2) (it was for this reason that the halves were cho- 
sen to overlap). Therefore, since \\TL\\ is bounded by 
VJ < rL d J, if we pick c in Eq. (1) sufficiently large, we 
will find that Tr[p(9, 0)He,o] - E a (0, 0) will also be small 
at 9 — 2ir; this follows from the statement that a small 
quantity multiplied by a fixed power of L is also small. 

Our claim, which we show in this section, is that 
Eq. (22) is satisfied by a p(9, 0) defined as follows for 
< 9 < 2tt. We pick 

dgp(9,0) = - J dT[A+( t r)-A-(- l r) 7 p(9,0)}, (24) 


where we define 



A+(it) = (27T)- 1 cxp[-(rA J B) 2 /(2<zlog(L))] x (25) 
dt d e H efi (t) c^[-(tAE) 2 /(2q log(L))]/(-*i + r), 



with q to be chosen later, and A~(— it) = (A + (ir))^. 
The time evolution of the operator dgTig t g'{t) is defined 
using the Hamiltonian Hg,-g, while the t dependence of 
A+- is defined via Eq. (25). 

To give some insight into the definition of A + , A~ , 
we note that if q were to be infinite, then they would 
project onto positive and negative energy parts of dgH 
at times ±ir, respectively. That is, for q = 00, we have 
A + {ir) = (2TT)- 1 J^dtdgHg.g^^/i-it + t). Let the 
matrix elements of the operator dgH in a basis of eigen- 
states of TC be written (dgH.) a b where the states have 
energies E a ,Et,. Let the states have energy difference 
— u = E a — Eb- Then, doing the integral over t we 
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find that A + (ir) has a matrix element between states 
a, b equal to exp[u)T](d$H) a b for — u> > and equal to 
zero for — w < 0. Similarly, A~(— it) has a matrix ele- 
ment equal to exp[— lot] {dgH) a b for —lo < and equal to 
zero for — lo > 0. Then, for any given time r, the inte- 
grand of Eq. (24) would be the same as that of Eq. (18) 
for 9 = 0, since in that case the only non- vanishing terms 
in Eq. (24) arc -A+(ir)p(0, 0) - p(0, 0)A"(-rr). 

What we will do later is to take a finite q instead. 
Physically, this means that rather than taking an adia- 
batic change in 9 which keeps us in the ground state, we 
instead "pass through" the level crossing when the gap 
closes at some 9^0, going from the ground state to some 
low energy excited state. 

Eq. (24) gives the change in p equal to the commutator 
of p with an anti-Hermitian operator, and hence gener- 
ates an infinitesimal unitary transformation of p. Thus, 
p continues to be a density matrix which projects onto a 
single state, defined to be ip(9, 0). 

As a first step, wc wish to show that for 9 = 9' = 
we can find a c such that d e p\(0, 0) - £\ i[S], p?(0, 0)] is 
small. We have d B (%(9, 0) - E 3 i[S$, P?(0, 0)] = 

d e p 1 (0,0)-d e p° 1 (e,-0)=Tr 2 [d ,p (p,e>)], (26) 

where all derivatives are evaluated at 9 = 9' = 0. To 
bound the right-hand side of Eq. (26), consider the trace 
of this term with any operator O with ||0|| = 1. This 
operator must be within half (1), so, using Eq. (18) to 
compute the derivative of p° with respect to 9', we obtain 
the expectation value Tr[Od0>p°(O, 9')] = 

oo 

d e >(0) = - j dT {Odg,Ho, e <(iT) + dg,H e ,e<(-iT)0). (27) 
o 

However, following the arguments from the previous sec- 
tion and the locality bounds, we can find a c such that 
Eq. (27) is small. In this case, the distance between O 
and dgiH is at least L/A — 2R, since de>H includes terms 
with x down to L/2 — R, while O is in half (1) so includes 
x up to L/4 + R. Note that {dg>He,e>) = at 6 = 6' = 0. 
Since we have bounded the trace of the right-hand side 
of Eq. (26) with all operators O with ||0|| = 1, we have 
bounded the trace norm of the right-hand side. 

B. Bound on Error Terms 

We now show that we can find a c such that the defi- 
nition (24) satisfies Eq. (23) in general. We wish to com- 
pute d e pi(0,O) - Pi (0,-0)] = e„(0) + e b {9) + 
e c (9). Here we define 

ciL 

e a = ~ J dTTr 2 [[A+(iT)-A-(-iT),p(e,0)}- (28) 
o 

[A+(iT)-A-(-iT),p°(e,-6)}]. 



ci L 

e b = - j Tr 2 p+(zr) - A~(-i T ), p°{6, -9)]] (29) 
o 

-d eP \{6,e'). 

e c = d e pl(9,9') - p?(0, -0)]. (30) 

3 

In Eq. (29), the derivative of p? is evaluated at 9 = —9'. 
We now consider each of these terms e a , eb, e c in turn. 

First, consider Eq. (28). In the definition of 
A + -~ as an integral over t, the integral over times 
\t\ > Ci(L/2 — R) has an operator norm bounded 
by \\deH\\f ltl>ci(L/2 _ R) eM-(tAE)y(2 q log(L)]. Thus, 
for any fixed q (to be chosen later) we can find a c such 
that this integral over times |t| > c\(L/2 — R) has small 
operator norm, and thus when commuted with p(0, 0) 
gives a term with small trace norm. 

Eq. (28) involves an integral of dgH(t) over time t in 
the definition of A + '~; we have shown that the con- 
tributions with times \t\ > c\(L/2 — R) may be ne- 
glected. Then, considering only contributions with \t\ < 
Ci(L/2 — R), we claim that, up to an error in the operator 
norm of order exp[— 0(L)], deH(t) can be written as an 
operator involving only terms not in half (2). That is, 
that \\deH(t) — Tr 2 [9eH(i)]|| is exponentially small in L. 
To show this, define U\2 to be the set of all sites j which 
lie in both half (1) and half (2); there are at most 2SV/L 
such sites. These sites are shown in the solid regions in 
Fig. (3). Define operators 0(t = 0) = 0'{t = 0) = d e H, 
and define the time evolution of O, O' by d t O = —i[0, H], 
while d t O' = —i J2^ Ul2 [O, Hi], i.e., the time evolution of 
O' includes only the sum over sites % which are either in 
half (1) or in half (2), but not in both halves. Then, using 
the arguments leading up to Eq. (4), we can show that 
for i e U12, the operator norm ||[0(t), Hi]\\ is bounded 
by 2J\\0\\g{cx{L/2 - R),L/2 - R), which is of order 
exp[— 0(L/2 — R)] for the given range of times t. Then, 
using the difference in the evolution equations for O, O' , 
we can bound \\0(t) - O' (t)\\ < £ ieC/l2 /* di'||[0, < 
2t{2SV/L)J\\0\\g{a(L/2- R),L/2- R). This quantity 
is also of order exp[— 0(L/2 — R)] for the given range of 
times t. Finally, we use the fact that Tr 2 [0'] = O' to get 
the desired result. 

From the above two paragraphs, it follows that 
up to small error in the trace norm, e a (9) = 

- °f dr [Tr 2 [A+ {it) - A~ (-zt)] , Tr 2 [p{9, 0) - p° {9, -9)]] . 
o 

Then, this is equal to the commutator of E\{9) with an 
anti-Hermitian operator. It generates an infinitesimal 
unitary rotation of E\{9) and therefore does not lead to 
any change in \E\{6)\. 

Next, consider Eq. (29). First consider the terms in 
the commutator involving A +, ~ acting on the left side of 
p°. As above, the operator dgH can be written in a basis 
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of eigenstates of 7i as (dgTL) a b, where the states have 
energies E a , E^. In dgHp°(9, —9) the only non-vanishing 
terms involve states with energy difference — w — E a — 
E\j > AE. Consider a matrix element (dgH.) a b with given 
u. This leads to a matrix element of A~ {— ir) equal to 
{deU)ab times 



v /2^glog(L) 
AE 



exp[-(TAE) 2 /(2q\og(L)} x 



00 di / 

— 9(u/) exp[-u/r] exp[- ? Iog(L)( W - W ') 2 /(2A£ 2 )], 

Z7T 



where we have converted the time integral to an integral 
in Fourier space. Since — u> > AE, Eq. (31) can be made 
small by choosing q sufficiently large. Thus, the trace 
norm of A~{— ir)p°(9, —0) is small, for all r > 0. Simi- 
larly, for A + , we find that we get a matrix element equal 
to {dgH) a b times 



v /2^1og(i) 



AE 



exp[-(TA£) 2 /(2glog(£)] x 



da/ 
~2tT 



6(-u/) exp[w'r] exp[- 9 log(L)(w - w') 2 /(2AE 



By choosing q sufficiently large, the integral (32) can be 
made equal to exp[wr], up to small error. Thus, the 
given matrix element can be made equal to {dgH) a b times 
exp[wr], up to small error. Therefore, the trace norm 

of - C /drTr 2 {[A+(iT) - d e H{iT)}p°{6, -6)} is small. 


These statements amount to saying that, with small er- 
ror in the operator norm, A + (ir) indeed is equal to the 
positive energy part of d$H(ir), while A~ is equal to the 
negative energy part. 

Now consider A + '~ acting to the right side of p°(0, —0), 
so that we consider p (0, —9)dgH. In that case, the only 
non- vanishing terms in (dgTi.) a b involve — u) = E a — Ef, < 
—AE. Repeating the argument above, we find that the 
trace norm of p°(9, — 9)A + {ir) is small, as is the trace 

norm of °f dTTr 2 {p (6, -0)[A~ (-ir) - dgH(-iT)}}. 


Therefore, up to small error, Eq. (29) is equal 

to - C fdT Tr 2 [d e n(iT)p°(9, -6) + p°(6, -9)dgH(-iT)] - 


dep\{6,ff), which equals - / dTTr 2 [<9 e W(rr)p°(0, -6>) + 



00 

p n {9,-9)d e H{-iT)] + JT dTTr 2 [d e H(iT)p°(9,-e) + 



p°(9,—9)dgH(—iT)}. This difference is equal to an inte- 
gral over t > c\L. For sufficiently big c, the trace norm 
of this integral can be bounded by any desired negative 
power of L. Thus, e& has small trace norm. 
Finally, consider Eq. (30). This is equal to 

dgp\{e,e')-dep\{e,-9) = T! 2 [d e ,p {e,e')], (33) 



where the derivatives are evaluated at 9' = —9. The trace 
norm of the right-hand side of Eq. (33) can be bounded 
by a negative power of L using the same arguments near 
Eq. (27), by considering an operator O that is entirely 
within half (1). The only difference to the arguments 
near Eq. (27) is that we compute the derivatives and 
expectation values at 9 = — 9' , rather than at 9 = 9' = 0. 
.gThus, using Eqs. (28,29,30) and picking sufficiently 
large q we find that Vg\E\{9)\ is small for sufficiently 
large c. A similar sequence of arguments permits one to 
bound Vg \E 2 {9)\. In the next section, we consider the 
expectation value of the translation operator on p(2ir, 0). 



VI. TRANSLATION OPERATOR 

Consider the operator tr((l, y){2, y)...(L, y)), which 
translates the sites with given y. The translation op- 
erator T which translates the entire system by one unit 
e^^l-jis the product of these loop operators over all y (there 
are an odd number of such loop operators). The ground 
state of Hofl is an eigenstate of T. If the ground state is 
non-degenerate, then it has eigenvalue ±1; without loss 
of generality we will assume in this section that is has 
eigenvalue +1. 

In this section we will show that the expectation value 
of T for p(2tt, 0) is opposite to that for p(0, 0), up to small 
error. We note that if Tig^o were to have a gap for all 
< 9 < 2n, then the results in this section would provide 
the last step in the topological argument discussed above. 
Instead, the results in this section will complete the ar- 
gument started in the previous section: p{2ir, 0) gives us 
a density matrix such that Trfp^Tr, O)?^^] — i?o(0,0) is 
small, but which, up to small error, has the opposite ex- 
pectation value for T . Since the difference in the expecta- 
tion of 7i is small, we can find a c such that the difference 
in expectation value decays faster than 1/L, and then we 
can find an Lq such that for L > Lq the state iP(2tt, 0) has 
an energy expectation value which is less than cln(L) /L 
above the ground state. However, since the expectation 
value of T is opposite for p(2ir, 0) compared to p(0, 0), up 
to small error, this state has an overlap on the ground 
state which is small. Thus, we will show in this section a 
contradiction under the assumption that the system had 
a gap AE which violated Eq. (1) and under the assump- 
tion that the system was translation symmetric, so that 
the ground state was an eigenstate of T. 

We first define a twisted translation operator, Tg.gi = 
Uy tr(exp[*§ a z }(l, y){2, y)...(L/2, y) exp[t^<j z ](L/2 + 
l,y)...(L,y)). Then, Tg t g> is a unitary operator and 
a symmetry of Hg^g/. Finally, given that T-0o(O,O) = 
V>o(0, 0), we have T e _ e ip (6, -9) = 4> {9, -9) for all 9. 

We will then show that dg\Tg fi ip(6,0) - tp(9,0)\ is 
small for all < 9 < 2tt. It will then follow that, up 
to a small error, Tr[p(0, 0)T ,o] = Tr[p(27r, 0)T 27r , ] = 



ft) 



— Tr[/9(27r, 0)To.o] = 1, thus showing that the twisted 
state indeed has the opposite expectation value for T. 
Here we have used the fact that for systems of odd width, 

T^nfi = — ^0,0- 

Consider the derivative d e (T efi ip(9,0) - ^(0,0)) = 

»^(i/^,0)- (34) 

y 

J dr [T efi A + {iT)-A+{iT) 
o 

-T efi A- (-it) + A- {-ir)]^(0, 0). 
This is equal to 

ci L 

{*E 5 (i,y) " / dr[To fi A + (ir)T^-A+(ir) (35) 

y o 

-T e , o A-(-iT)T d -i+A-(-iT)}}iP(e,0) + 
^J2 S a-v)- [ dr[Te,oA + (ir)T e ^ 

y o 

-Te,oA-(-iT)T d -i}}{T e ,oi>(6,0) - ^(0,0)} 

The last term of Eq. (35) is equal to an anti-Hcrmitian 
operator acting on Tg yO tp(9,0) — tp(9,0), and thus does 
not change the norm of this state. Thus, we need to 
bound the norm of the first term. This term is equal 
to an anti-Hcrmitian operator, that we define to be Oi, 
acting on tp. The norm square of this term is equal to 
Tr[O\Oip(0, 0)]. As shown in the previous section, up 
to small error in the operator norm, A + '~(ir) can be 
written entirely as operators in half (1). Therefore, 0\ 
can be written entirely as an operator in half (1); that 
is, the operator norm \ \0\ — Tr 2 [Oi]|| is small. Thus, the 
norm square is, up to small error, Tr[O|(9ipi(0,O], which, 
again up to small error, is equal to Tr[0\Oip ( l(6, — 0)], 
since \ Pl (0,0) - p\ (0,-0) | is small. 

We claim, however, that this last expectation value 
is small. To show this, consider dg(Tg t -gtpo(9, — 0) — 
ipo(9,—9)). This is equal to zero. However, this 
derivative can be written as an operator O act- 
ing on Vo, with O = iJ2y( S (i, y ) - S {L/2+i, v )) - 
r o °° dT(T e _edgHe,e'(iT)T^ e - d e TLg y g-{iT) 

Tg-gdgiHefi'{iT)Tg]_ g + dg'Hg ft> (it)) , where the deriva- 
tives are taken at = —0'. Since Oipo = 0, it follows 
that Tr[Ot<V(0,-0)] = 0. However, up to small error, 
O = 0\ + 02, with 0\ the operator considered above 
and O2 defined to be a similar operator acting only in half 
(2). Then, 2Tr[OjOip°(0, -0)]+2Tr[O|O 2( o°(0, -0)] = 0. 
However, using the locality bounds, the second term can 
be made small for large enough c [23], and thus the first 
term, Tr[O|Oip?(0, -0)], is small. Therefore, we have 
shown the desired result. 



VII. DISCUSSION 

The main result is Eq. (1), obtaining a bound on the 
energy gap for spin models in arbitrary dimensions. In 
order to obtain this result, we have introduced a set of 
loop operators, and proven a bound on connected cor- 
relation functions. This bound on correlation functions 
did not rely on the system being a spin system; rather, 
it is valid for any Hamiltonian such that the Hi have 
bounded operator norm, and such that the interaction is 
finite range. Below, we generalize this bound on correla- 
tion functions to certain other systems as well. 

We note that for the case of higher spin representations 
of SU(2), Eq. (1) follows automatically from the result 
for spin- 1/2, so long as the total spin within all unit cells 
is half-odd: the higher spins can be written as various 
combinations of spin-1/2 spins, and if the total spin in the 
unit cell is half-odd then there will result an odd number 
of spin-1/2 spins in each unit cell. Suppose, for example, 
a unit cell contains one spin-1 spin and one spin-1/2 spin, 
giving a total spin of 3/2 which is half-odd. Then, the 
spin-1 can be written as two spin-1/2 spins. Let these 
two spins be called S± , S2 and let the Hamiltonian in- 
clude only terms symmetric under interchange of Si, S 2 . 
This new Hamiltonian has three spin-1/2 spins in each 
unit cell, and hence falls within the class of Hamiltonians 
considered above. Then, there are two different sectors 
of the Hilbert space with no terms in the Hamiltonian 
coupling these two sectors: one sector in which Si,S2 
form a spin-0, and one in which they form a spin-1. By 
adding a term coupling Si to S 2 to the Hamiltonian with 
a sufficiently large, negative (ferromagnetic) coefficient, 
we can ensure that the ground and first excited states lie 
in the sector in which Si, S2 have total spin-1. Then, the 
existence of a low- lying state satisfying Eq. (1) for the 
new system with only spin-1/2 implies the existence of 
such a low-lying state for the original system with both 
spin-1 and spin-1/2. It would also be interesting to gen- 
eralize these results to other groups SU(N), as well as to 
consider the case of even V/L. 

We finish with two conjectures. First, we conjecture 
that the same Eq. (1) holds for systems with an even 
width, so long as the width V/L is of order L d ~ Y and so 
long as d > 1. For d = 1, this result is of course not true, 
as Haldane gap behavior is possible. 

Second, consider the thermal expectation value of 
T at an inverse temperature (3, defined by (T)p = 
Tr[exp(-/3W)T]/Tr[exp(-/3W)]. We conjecture that 
there is a constant c, depending on TL, d such that for 
(3 > c\og(L)/L the given thermal expectation value {T)p 
vanishes in the limit L — > 00 for systems of odd width. 
We base this conjecture on the following physical obser- 
vations: for ferromagnetic systems, there are spin wave 
excitations, with dispersion relation E oc k 2 . It may 
be shown that the presence of these excitations causes 
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(T)p to vanish for (3 of order L as L — > oo. For antiferro- 
magnetic systems, the translation symmetry is broken by 
the antiferromagnetic ordering (in fact, for these systems, 
the true ground state has translation symmetry and is a 
superposition of different broken symmetry states, but 
there are low-lying states with different expectation val- 
ues of T so that for (3 of order L, the expectation value 
(T)p vanishes). Finally, for spin liquid systems, there is 
a low-lying excited state with the opposite expectation 
value of T compared to the ground state, as we have 
found above. We leave a proof of both of these conjec- 
tures for future work. 



APPENDIX A: MARKOV PROCESSES AND 
LOCALITY 

Consider a system with a probability pi of being in 
state i and a transition matrix Ty so that the equation 
of motion is dtPi — ^ijPj ■ For the total probability 
to be conserved, we have J2i Tij = 0, which guarantees 
that Tij has at least one zero eigenvalue. Let us assume, 
further, that all eigenvalues of are real. This includes 
all systems for which the stationary state (given by the 
zero eigenvector of T^ ) obeys detailed balance. A typical 
example of such a process would be the Monte Carlo 
dynamics of a statistical mechanics system. We will first 
derive a suitable generalization of the locality result (13) 
to systems governed by such a Markov process, and then 
discuss the implication for statistical mechanics systems. 

Let us assume that the spectrum of T is such that there 
is only one zero eigenvalue, with right eigenvector and 
that all other eigenvalues A are negative with A < —A, 
for some A > 0. Assume p { \ is normalized by J2iPi = 1- 

Then, introduce various quantities to be measured, 
A,B, so that the expectation value of A is given by 
(A) = J2i^iPi- We can write this slightly differently 
by introducing for each such quantity a diagonal matrix 
given by An = Ai for all i and Aij = for i j. Fur- 
ther, introduce an additional vector Jj, such that 7j = 1 
for all i. This vector 7j is a left eigenvector of T with zero 
eigenvalue, since J2i Tij — as mentioned above. Then, 

(A) = 

We can now consider expectation values of quanti- 
ties at different times: (A(t)B(0)) = PAexp[Tt]Bp° = 
P enp[-Tt}Aexp[Tt}Bp = PA(t)B(0)p°. In these 
equations, P denotes the transpose of the vector I (I 
is real, so no complex conjugation is necessary), and we 
have left off all the indices on vectors I,p and matrices 
A,B,exp[±Ti\: the product is evaluated following the 
usual rules of matrix multiplication. In the sequence of 
equalities above, the first equality defines the time evo- 
lution of the system, the second equality follows since 
Jexp[— Tt] = I, and the last equality follows since we de- 
fine A(t) by the equation of motion: d t A(t) = [A(t),T]. 



It is then possible to extend this definition to operators 
separated by an imaginary time separation: (A(ii)B(0)). 

Now, consider a typical physical example: an Ising sys- 
tem, governed by Monte Carlo spin flip dynamics, with 
A(0),B(0) representing the value of two different spins 
which are separated in space. In such a case (as well as 
in many others), it is possible to obtain a bound similar 
to Eq. (4). Assume that the matrix T can be written as 
a sum of matrices Tj, with finite interaction range R and 
with a bound ||Tj|| < J, where i is a site index. Define 
f(t) = ([A(it), -B(O)]). Since f(t) < \\[A(it),B( 



fit) < 



2N B S\\A\\\\B\\\2JSt\ l / 2R e 2JS W 



{l/2R)\ 
= N B 



(Al) 



B\\g(t,l). 



At this point, from the existence of a A > and a spec- 
tral representation with all eigenvalues real [24] follows a 
result similar to Eq. (13): 



\(A(0)B(Q)) c \ < i-2JV B 
2(1 



B\\g( Cl l,l)+ (A2) 

B n- C1 AEl/2_ 



Therefore, if there is a Markov dynamics that gives rise 
to the equilibrium probability distribution p® which has 
a A > 0, then there is an exponential decay of correlation 
functions in space. An example is a spin system in the 
paramagnetic phase with Monte Carlo spin flip dynam- 
ics. The converse is not necessarily true: a spin system 
in the paramagnetic phase with spin exchange dynamics 
does not have a A > but instead has spin correla- 
tions which decay with a power law time. However, this 
dynamics gives rise to the same equilibrium probability 
distribution as the spin flip dynamics does, and hence 
has exponentially decaying correlations in space. 
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FIG. 1. Illustration of the bound on |/+(0) - / + (0)], as 
described in the text. The vertical line describes a J-function 
spike in f(u). This produces a Gaussian in f(ui). The integral 
over u> of the Gaussian is the same as the integral of the 
5-function; however, the shaded region of the curve of the 
Gaussian falls above lu — 0, and hence does not contribute to 
/+(0). This leads to the difference between /+(0) and / + (0). 
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x=0 x=L 

FIG. 2. Plot of the system, showing the ^-coordinate along 
the length axis. The x is shown ranging from x = to x = L; 
due to the periodicity of the system, x = is identified with 
x = L. The y coordinate specifies the position in the direc- 
tions normal to the length, as well as specifying the particular 
site in each unit cell. The twist angles are noted; the twist 9 
changes the boundary condition near x — L, while the twist 
9' changes the coupling between sites near x = L/2. 
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x=0 x=L 

FIG. 3. Plot of the system, showing the twists and coordi- 
nates as before. The halves of the system have been shaded 
in. The shading at the left and right side of the system (di- 
agonal lines going up and right) denotes sites in half (1), the 
shading in the middle (diagonal lines going up and left) de- 
notes sites in half (2). The solid shading denotes sites in both 
halves; the length of the solid region is at least 27?, so that 
the Hamiltonian can be written as a sum of terms, each of 
which is contained in only one half. 
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